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Thermoelectric properties of disordered systems 
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1. Introduction 

The electronic properties of disordered systems have 
been the subject of intense study for several decades. 
Thermoelectric properties, such as thermopower and 
thermal conductivity, have been relatively neglected. A 
long standing problem is represented by the sign of the 
thermoelectric power. In crystalline semiconductors 
this is related to the sign of the majority carriers, but 
in non-crystalline systems it is commonly observed to 
change sign at low temperatures. In spite of its appar- 
ent universality this change has been interpreted in a 
variety of ways in different systems. We have developed 
a Green's function recursion algorithm^"' based on the 
Chester-Thellung-Kubo-Greenwood formula^"'*-' for cal- 
culating the kinetic coefficients Ly on long strips or bars. 
From these we can deduce the electrical conductivity cr, 
the Seebeck and Peltier coefficients 5 & 11 and the ther- 
mal conductivity k, as well as the Lorenz number Lq- We 
present initial results for ID systems. In ID we observe 
a Lorentzian-like distribution for the thermopower which 
is modified by the presence of inelastic scattering. This 
could give rise to non-negligible quantum fluctuations in 
macroscopic systems at low temperatures. 

Within the linear response approach, the responses of 
a system to an external electric field E and a temperature 
gradient VT up to linear order^^ are 



(j) = |eri(|e|LnE-Li2T-iVT) 

= crE - aSVT 

(j,) = e-2(|e|L2iE-L22T-VT) 

= STO) - KVT. 



(1) 



Here, (j) and (j^) denote charge and heat currents, 
i,j = 1,2 are the kinetic coefficients and e, T, 
(T, S, and K are electron charge, temperature, conduc- 
tivity, thermopower and heat conductivity, respectively. 
Other commonly considered thermoelectric transport co- 
efficients are the Peltier coefficient, 11, and the Lorenz 
number, Lq 
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with /cb the Boltzmann constant. Calculations of the 
thermoelectric transport properties that have been done 
previously are based either on the CTKG formulation^"''^ 
or other perturbative methods.^' These methods are 
macroscopic in nature and rely on the explicit, phe- 
nomenological form a{Ep,T). A microscopic approach 
to determine the thermoelectric transport properties 
without assuming explicitly a functional form of a would 
be quite instructive. To our knowledge, there has only 
been Ref.^^ which reported on the distribution of 5* in a 
ID disordered wire in the localized regime. The authors 
were able to analytically and numerically show that in 
this case the distribution of S is Lorentzian. So far there 
has been no other microscopic calculation of S, 11, K or 
Lq in disordered systems near the MIT. 

2. The recursive Green function method 

This method was used previously to calculate the dc 
and ac conductivity tensors and the density of states of a 
d-dimensional disordered system^' with and without a 
magnetic field. Here we will show briefly how the method 
can be adopted to calculate all kinetic coefficients and 
thus the other transport properties. 

The ID Hamiltonian can be written as 



TV 



H 



(5) 



i=i 



where the disorder terms tj G [— W/2,iy/2] are ran- 
domly distributed. In adding one more site to the chain, 
the Hamiltonian can then be considered as an unper- 
turbed Hamiltonian of the first — 1 sites and the de- 
coupled A^*'' site plus a perturbation.^^ Equation (5) 
then simplifies to 



-iV-l 



(6) 



The last part in the Hamiltonian, iJ', is the coupling be- 
tween the (A^ — 1)*'' site and the A^*'' site. This coupling 
can be treated iteratively using the one-particle Green 
function G'^(z^) as 



(7) 
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where = E ± i^, where E is the energy and ^ may 
be interpreted as an inverse hfetime, e.g. due to inelastic 
scattering. 

3. Results and Conclusions 

For brevity, we shall restrict ourselves here to studying 
the behavior of 5 (^12). In complete analogy to the 
formulation for a (= Lii),^'^°) we can compute Li2{T = 
0) using the recursive equations 



= {E-n)& 



(JV) 



= 9n^ [{XN-XN-lf9M-l + ^^^-l] 

ft, = ft-i-'^ati [{xN-XN-i)'^gN-i ^ 



(8) 
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The initial condition is fo=hQ= 0. From Eq. (3) we 
see that 



|e| 1^' 4 Siv J 
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(9) 



where sjv = ttKN aa'^^^ / Ae^ . The distribution of L12 (in 
units of e^/wha) is shown in Fig. 1. 
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compared with the results of the recursive approach as 
shown in Fig. 2. 
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Fig. 2. Distribution function P{t) for zero and finite value of ^ 
and the same parameters as Fig. 1. The dashed lines correspond 
to a Lorentzian (left) and a Gaussian (right) fit. 



Using 0{N) arithmetic operations, where is the 
number of sites in the system, we can compute itera- 
tivcly all kinetic transport coefficients Lij. The maxi- 
mum value of N is mainly limited by the computer time 
available and the intended precision of the results. In 
this short note, wc have considered only the transport 
properties of a chain initially having N sites of a ID 
tight-binding model with nearest neighbor hopping in 
the limit T = 0. It is straightforward to extend the 
computations to higher dimensions^^^ and results will 
be presented elsewhere. Detailed derivations of relevant 
equations and recursion relations arc given elsewhere. -'^^^ 

We have shown that the distribution of the dimension- 
less thermopower r is indeed Lorentzian as predicted^^ 
but that this changes to Gaussian on introduction of 
a finite lifetime. Similarly the distribution of L12 is 
Lorentzian-like (albeit asymmetric) for ^ — > but also 
becomes Gaussian for finite ^. 



Fig. 1. Distribution function P(Li2) for zero and finite value of ^ 
at E = jj, = —1, W = 2 and T = 0. Data is computed from 
10000 samples and length Na = 3 x 10® a. P is normalized to 
1 at the maximum. The dashed line correspond to a Gaussian 
(right) fit. 

A further interesting quantity is the dimensionless 
thermopower t defined via 



dg{E) 



2Trg{E) dE 



E=Et: 



2^dE 



In 



(10) 

with A the mean level spacing near Ep and g the con- 
ductance. In the limit of weak scattering the domi- 
nant effect of the disorder is simply to introduce a finite 
mean free path at energies well inside the band. This 
implies that on the average (10) can be approximated 
by r = -fTr [ReG+{E)]. The full distribution can be 
shown to have a simple Lorentzian form^^ which can be 
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